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Abstract. Painleve test of the coupled Gross-Pitaevskii equations has been 
carried out with the result that the coupled equations pass the P-test only 
if a special relation containing system parameters (masses, scattering lengths) 
is satisfied. Computer algebra is applied to evaluate J = 4 compatibility 
condition for admissible external potentials. Appearance of an arbitrary real 
potential embedded in the external potentials is shown to be the consequence 
of the coupling. Connection with recent experiments related to stability of two- 
component Bose-Einstein condensates of Rb atoms is discussed. 



PACS numbers: 03.75.Fi,05.45,32.80.P 
1. Introduction 

Recently there has been a growing interest in the Gross-Pitaevskii (GP) equations 
EI describing two-component Bose-Einstein condensates (BEC) in external trap 
potentials [3-22]. In the absence of the confining potential, the GP equations reduce 
to the coupled non-linear Schrodinger (NLS) equations which play an important role in 
optics Coupled GP equations are also used to describe Josephson-type oscillations 
between two coupled BEC [9-11], spin-mixing dynamics of spinor BEC [12-15], or 
to explore such interesting field of matter waves as possible atomic soliton lasers 

An efficient tool of the analysis of the non-linear partial differential equations is 
the Painleve (P) method TP which serves to explore the singularity structure of 
the underlying equations, and establish integrability conditions |28j. The P-analysis 
of the single NLS equation has been performed by Steeb et al [221, and the damped 
NLS (or the GP) equation has been investigated by Clarkson [20]. A fairly large 
class of coupled NLS equations including third order dispersions have been analyzed 
by Radhakrishnan et al |3T]- Recently the symmetrically coupled higher-order NLS 
equations have been tested by using the P-method |32| . 

Because of the experimental developments in forming two-component BEC |16) 
and the possibility to confine BEC in a linear shape we shall perform the P-test 
of the coupled one-dimensional GP equations in order to establish certain necessary 
conditions of integrability. (The term integrability is used here in the general sense 
| 27l 121] involving P-property and soliton formation.) The results obtained for the 
trap potentials are similar to those found by Clarkson [301 in the case of the damped 
NLS equations: the trap potential should be linear and/or quadratic in the coordinate 
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variable x. In the quadratic case a source term depending only on time t should also 
be present in the external potential V{x,t). 

A novel feature of our analysis is the possibility of the appearance of an arbitrary 
common potential term V{x,t) within the confining potentials Vi{x,t) and V2{x,t). 
Its presence may prove useful for fine tuning experiments with two component BEC. 
We consider the system of coupled GP equations in its most general form containing 
different masses, external potentials, and mutual coupling strengths. As a result 
we shall derive compatibility conditions, the fulfillment of which depends on the 
parameters characterizing the GP equations. We show that in a particular experiment 
U] employing two hyperfine states of Rb atoms as components of BEC, the vortex 
stability corresponds to just the parameter ratios satisfying our general formula derived 
in this paper. 

The organization of this paper is as follows. In section 2 the P-analysis of two 
coupled GP equations will be carried out including the determination of the leading 
orders, the recursion relations, the resonances, and the compatibility conditions. 
The consequences of the compatibility relations for the potentials are discussed in 
section 3 where also other consistency requirements are studied. In section 4 we 
make comparisons with earlier results and investigate compatibilities with existing 
experimental and numerical findings related with two-component BEC. Section 5 is 
devoted to a short summary. 



2. Painleve test 



Let us consider the following (1 + 1) dimensional inhomogeneous NLS equations for 
the wave functions V'l: '4'2 with the external potentials Ui{x,t), U2{x,t) 

\h—i;i{x,t)={ + Ui{x,t) + Uii + Ui2 |V'2(x,t)|' Mx,t) + C/io, 

at \ 2mi / 

(2.1a) 

i;i— V'2(a;,t)=f - — + C/2(x,t) + U21 \Mx,t)f + U22 |V2(x,t)n ^2(2;, i) + U20 

(2.16) 

which, in the absence of the inhomogeneities Uiq and C/201 £^re commonly called the 
coupled Gross-Pitaevskii equations PEl- 

Here rrii denotes the mass of the atomic species i (i = 1, 2) of the two-component 
BEC gas and Uij is related with the interactions between the atoms i and j (i, j — 1,2) 
via the relation Uij = 2'Kh^aijNj / A where Nj means the number of atoms in the 
jth component of the BEC, is the scattering length characterizing the interaction 
between atoms i and j, A represents a general cross sectional area confining species i 
and j, and = rriimj/ {mi -t- rrij) is the reduced mass. 

By introducing the new parameters 

A = ^, T., = if/„ (,,.1,2), (2.2a) 

and notations 

u = ^i, w = ij2, V^ = lu,, V,o = lu,o, (i = l,2) (2.26) 
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we write the GP equations into the standard form of the P-analysis 

iuf + Xuxx — Til u — Ti2 \w\'^ u — Viu + Vio, (2.3a) 

iwt + 'dw^x - T21 \u\'^w - T22 |wp w= V2W + V20 (2.36) 

where Tij and A, ■& stand, as defined by equation (|2.2a|l . for the interaction and mass 
parameters, respectively. 

In order to apply the P-analysis, we first complexify all variables to obtain 
equations ^2.'Sai -b) in the form [v = u*, z = w*): 

iut + Xuxx — Tiiu'^v — T12WZU =Viu + Vio, (2.4a) 

-ivt + Xvxx - Tiiuv"^ - T12WZV = ViV + ViQ, (2-4&) 

iwt + -dwxx - T21UVW - T22w'^z — V2W + V20, (2-4c) 

-izt + -Dzxx - T21UVZ - T22wz^ = V2Z + 1^2*0 (2.4rf) 

where the functions u, f , w, z are treated as independent complex functions of the 
complex variables x and t, and Vi {x,t), Viq {x, t) , V2 {x, t) , V2Q {x, t) are formal complex 
conjugates olVi{x,t), Vio(a;,t), V2{x,t), V2o{x,t), respectively. 
The next step is to seek the solutions of (I2.4al -d') in the form 



00 

u{x,t) = (f{x,t)^u.j{t)(j)^{x,t), v{x,t) = (l)'^{x,t)^Vj{t)(l)^{x,t), (2.5a) 

3^0 j-O 



00 00 



w{x,t) = (j}''{x,t)2_^Wj{t)(t)'{x,t), z{x,t) = (t)''{x,t)2_^Zj{t)(j)^x,t), (2.5&) 

with the Kruskal ansatz 

cj){x,t)^x^^{t), (2.6) 

and ^(t), Uj{t), Vj{t), Wj(t), Zj{t), j = 0, 1,2, . . . being analytic functions of t in the 
neighbourhood of a noncharacteristic movable singularity manifold defined by = 0. 
Similarly, the external potentials Vi confining specimen i is also expanded about the 
singularity manifold = as follows {i = 1,2) 

V{x,t)=f2V,,it)cl^^ix,t), v,At)^U^^^^^) . (2.7) 

Substituting expansions ^2.5a\ -b) and (|2.7fl into equations ^2Ad -d) and equating 
like powers of (f) we obtain: 

i. ) equations for determining the leading orders p, g, r, s; 

ii. ) recursion relations for deriving the functions Uj, Vj, Wj, Zj. 

In order that equations (|2.4al -rf') pass the Painleve test it is required that the numbers 
p, q, r, s be non-positive integers. Moreover, the recursion relations should be 
consistent in all order of j including the resonances. 

2.1. Determination of the leading orders 

To determine the leading orders p, q, r, s appearing in the expansions (|2.5al -6'). it is 
sufficient to consider the expansion upto the zeroth order, j = 0. By substituting this 
truncated version of expansions ^2.b(A -h) into H2.4af) we obtain 

mo,tr + mopr~^cl)t + Xuopip - - Tiiulv^cf^'P+'i - Ti2U^woZacj)P+-'+' 

= Fiuo^^' + Vio. (2.8a) 
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Three completely similar expressions arise from the substitution of the truncated 
version of expansions (|2.5al -&) into the remaining three equations Ij2.46l -(i): 

- ifo,t0'' - iwo#''"Vt + Avo9(9 - 1)0''"^ - ThWqUo^^'^^^ - Ti2VoWozo(t>'^'^''~^'' 

= V{vo<l>'' + V{o, (2.8&) 

iwo,t<f>'' + ™o^'/'''"Vt + ^wor{r - 1^-^ - T2iUoVoWor^''+'' - T22wgzo0"'+^ 

= ^2^00" + V20 , (2.8c) 

-izo,t^' - izoS0^-Vt + ^Z0S{S - 1)0^-2 _ T2iUoVoZor+''+' - T22t«0^0 

= V2*zoq^' + V2I. (2.8rf) 

By demanding the leading order terms of equations ^2.8a\ -d) to vanish one obtains 
the following equations 

Xp{p - 1) = TiiUqVq + T12W0Z0, 

Xqiq - 1) = TiiUqVq + TuWqZq, 
■dr{r - 1) = T21U0V0 + T22W0Z0, 

Z?s(s - 1) = 121^0^0 + T22W0Z0 

and 



P + q=-2, 
r + s = -2 

from which the leading orders can uniquely be determined to be 

p ^ q = r = s = —1. 
For later use we infer from equations ^'2.9d -d) the useful relation 



(2.9a) 
(2.96) 
(2.9c) 
(2.9rf) 

(2.10a) 
(2.10&) 



UqVq 
WqZq 



(2.11) 
(2.12) 



2_ f T22 -Ti2\ fx 
A {~T2i Tn J {1} 

with A = T11T22 — Ti2r2i. If accidentally A = happens then we may use the relation 
UqVq/woZq = const instead of H2.12|l . which case needs a special consideration. 



2.2. Recursion relations 

The next step of the P-analysis is to again substitute expansions ^2.ba\ -b) and H2.7|l 

with the leading orders p = q = r = s = —1 into equations ]2.4:aH -d). After some 
algebra we obtain the recursion relations 

/ Ql -Tiiul ~Ti2UqZo ~Ti2V 

-Tiiv^ Ql -T12V0ZQ -T12V0W0 

-T21VQWQ -T21UQWQ Q2 -T22WI 

y-T2iWo2:o -T'2iUoZo -722^0 Q2 



\ 


1 




FA 












Wj 













(2.13a) 



where j = 1, 2, . . . and 

Ql = X{3 - - 2) - 2TnMoVo - ^12^0^0, (2.136) 
Q2 = i?(j - l)(j - 2) - 2T22WQZQ - T21U0V0, (2.13c) 

i-i 

F-j = -iWj-2,t - i(i - 2)Uj^i(f)t + ^ {TiiUmUj-rnVo + Ti2U™WoZj-m) 

m— 1 

+ ^ ^ {TiiUmUl-rnVj-l + Ti2UmWj_(Z;_„) + ^ Vl,iMj-;-2 + • (2.13d) 

i=l m=0 i=0 
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Here we use the notation that whenever an index is less than zero, then the expression 
itself is zero (for example, Vioj_3 = for j < 2). Furthermore Gj is obtained from 
Fj by interchanging ui and vi and letting i — i, Vi^i — > V^i and Vlo,; ~^ The 
expressions Hj and Kj can be obtained, respectively, from Fj and Gj by interchanging 
ui and w;, vi and z/, T12 and r2i, Tii and T22 , and letting Vl ^ ^2, Vw ¥20- 

The expressions Fj, Gj, Hj, Kj at a given j depend only on the expansion 
coefficients ui, vi, wi, zi with I < j. Therefore the equation H2.13q|I represents recursion 
relations for the determination of the unknowns Uj , Vj , Wj , Zj from the knowledge of 
the prior calculated coefficient functions ui, vi, wi, zi with I < j. 



2.3. Resonances 

The above recursion relations H2.13Qjl determine the unknown expansion coefficients 
uniquely unless the determinant of the matrix Q{j) is zero. Those values of j at which 
the determinant det{Q{j)) becomes zero are called resonances. After some calculation 
one obtains 

det(g(,)) = XHHj + 1) fU - mj - 4) - 3, + 4 - ^ '^T^^-oVo + m.woZo ^^ 

(2.14) 

so that the resonances of the coupled GP equations H2.3al -&) are as follows 

j,es = {-1, 0,0,3,3, 4,ji,j2}. (2.15) 

Here ji and j2 are the roots of the expression contained in the last parentheses of 
equation (|2.14() and can formally be given as 



3 1 / 'dTiiUpVo + XT22W0Z0 „ 
Ji,2=2±2V8 ^ 7 eZ. (2.16) 

As indicated, the resonances ji and j2 must be integers so that the square-root 
should be odd integers. From this one gets a condition 



dTiiUgVo + XT22WQZ0 

— 7=2m+l, m = 0,1,2,... (2.17) 

involving the characteristic parameters Tij, A, d of the GP equations ^2.'Aa\ -b). The 
number m can be considered as a classification number which classifies possible 
external potential families for which the system l^'2.3a\ -b) is integrable (in the general 
sense of integrability TP, ) . 

By re-arranging H2.17|l and using relation (|2.12(l one gets a more explicit condition 
necessary for any coupled GP equations to pass the P-test: 

2TnT22-(./A)TnT,2-(A/^)T2,r22 ^ J_ ^ . 0, 1, 2, . . . (2.18) 

T11T22-T12T21 16 ' J' . > ' V y 

It is also clearly seen that this expression depends only on the ratios X/d, T11/T21, 
and T12/T22 involving the characteristic parameters of the GP equations. 

In summary, any coupled system of GP equations H2.3al -61 passes the Painleve 
test only if its characteristic parameters A, t9, Tij — 1, 2) obey the relation (|2.18|) . 
otherwise it is probably not integrable. (See discussions about connection of P-test 
with integrability in references |22ii2Sj.) 
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2.4- Compatibility conditions 

At each element of jrcsi the recursion relations H2.13a|l cannot be used for the 
calculation of the expansion coefficients. At these indices arbitrary functions may 
arise in the expansions l|2.5al -6). However, in order that the solution be expressible 
in the form of the expansions H2.5qI -M and (|2.7(l . the recursion relations should be 
identically satisfied at j € jres ■ The investigation of these specific requirements leads 
to relations called compatibility conditions which impose restrictions for the external 
potentials Vi{x,t) i — 1,2. We note that only the positive resonances are of interest. 



2.4- 1- Compatibility condition belonging to resonance j = 3. Let us consider 
equations p.lHajl at j — 3 and use equation l|2.12|l . The result is an equation 

-Tiiul -T12U0Z0 -T12U 

-TuUoVq -T12V0Z0 -'T12V0W0 

-T21U0W0 -T22W020 

'T21U0Z0 -T22Z0 -~T22WoZi: 



f-TiiUoVo 
-T21VQW0 

















G3 




W3 






/ 


\Z3J 







(2.19) 



0(3) 



whose matrix 0(3) has rank two. Indeed, by multiplying the first row with vq, the 
second with uq, one gets a matrix which possesses identical elements in its first two 
rows. Performing similar manipulations, one can make the third and fourth rows also 
identical. It then follows that the above recursion relation can only be consistent if 
the following compatibility conditions hold 



= G3-U0, 
H3Z0 = K3W0 . 



(2.20a) 
(2.206) 



We emphasize that the above conditions are not independent from each other because, 
for example as shown by l|2.13d|l . F3 contains elements Wi, Zi with i < 3. Similarly, 
G3, H3 and K3 also contain all types of expansion coefficients Ui, Vi, Wi, Zi with i < 3. 



2.4.2. Compatibility condition belonging to resonance j — 4. By taking the recursion 
relations (|2.13ajl at j = 4 and applying equation H2.12|l . one arrives at the following 
equation 

-T12U0Z0 

-T12V0Z0 

M - T22W0Z0 
2 



/4A 



T21V()Wq 
V -T21V0ZQ 



-Tiiul 
4A - TiiUqVo 
-T21U0W0 
—T21UQZQ 



-T12U0W0 
-T12V0W0 
-T22WI 



-T22Z0 -T22W0ZQJ 



\ 


/uA 




(FA 








Gi 




Wi 




Hi 


J 


\zi) 




\Kil 



.(2.21) 



0(4) 

In the general case the matrix 0(4) has rank three which means that only three of its 
rows are independent. Using this fact, after some calculation we obtain the following 
compatibility condition 



T2 



L [PiVfi 



G4U0 



> [lliZQ 



KiWQ 



0. 



(2.22) 



We should investigate also the possibility when rank((3(4)) = 2. In this case the 
compatibility condition decomposes into two distinct parts as it can be seen in the 
following way. The rank of a matrix equals the maximal order of its non-singular 
submatrices. We should thus calculate the determinants of all third order submatrices 
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of Q(4) and investigate the cases when they simultaneously become zero. After a 
simple but lengthy calculation the following results are obtained for the determinants 
of the four third-order submatrices of the matrix (5(4): 

|l6At?T2iUozo ; — 16Ai?ri2WoWo ; 16Ai9Ti2Wo-2o ; — 16A'!9T2iMq| . (2.23) 

Because A and are the non-zero mass parameters [sec definitions (|2.2a|l ]. it is clear 
that the subdeterminants vanish only if T12 = T21 = 0. This situation however 
corresponds to two uncoupled GP equations. We thus arrived at the compatibility 
condition found by Clarkson j^O] for the single GP equations: 

FiVo + G4U0 = 0, (2.24a) 

H4Z0 + K4W0 = . (2.246) 

We note however that the general compatibility condition, as given by equation 
H2.22|l . leads to more complicated external potentials (discussed in the next section) 
than that obtainable from equations l^'2.24a\ -b) with T12 = T21 = 0. 

2.4-3. Compatibility condition belonging to resonances ji and j2. We now consider 
the matrix Q{ji^2) with ji^2 taken from equation (|2.16|) . In general the matrix Q{ji,2) 
has rank three from which the following compatibility condition arises 

■dwQZo(^Fj^ .2VQ + Gji 2^0^ - Awowo^-ffji.s^o + ■^^^1,2^0) = 0. (2.25) 

As before we analyze also the case when rank(Q(ji_2)) — 2. After a lengthy but 
simple calculation we obtain the determinants of the four third-order submatrices of 
the matrix Q{ji2) to be: 

|r2iA^ {uovof (wozof zq ; -T2iA^ (uovaf {wozof ; 

T2iA^ (uovo)^ {wqZo)'^ zq ; -Ti2A^ {uoVoWQZofY (2.26) 

Now, as clearly seen the determinants (I2.26f) vanish simultaneously only if 
T12 = T21 — (or A = which case is not considered here). On the other hand, 
for this decoupled case one can determine the values ji and j2 by using definition 
H2.16|l and relation H2.12|l to be ji = 4, j2 = —1. But then, as it can be checked easily 
by using (|2.13Qf) . the corresponding compatibility conditions reduce to those already 
discussed in connection with equations H2.24qI -6'). We note however that, depending 
on the experimental situations, it is possible to obtain resonance values ji and j2 
greater than four. We should then use equation ()2.25f) for drawing conclusions about 
the admissible form of the external potentials. 

3. Possible forms of the external potentials 

In the preceding section we have found equations, called compatibility conditions, 
that must be fulfilled in order that the GP equations pass the P-test. In this 
section we exploit the consequences of these equations for the general form of the 
external potentials Vi, Vio, V2, V20 appearing in equations (|2.1al -6') and H2.3al -b). The 
experimental realization of such potentials may lead to detection of stable structures 
(like vortices) in BEG. 

Although the compatibility conditions are related with indices j at which the 
recursion relations ()2.13q|) do not apply to the calculation of the unknown coefficients 
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Uj{t), Vj{t), Wj(t), Zj(t), the equations (|2.2UqI -6). (|2.22|) and (|2.25|) can be reduced 
to expressions in which only the zeroth order coefBcient functions uo(i), ^^{t), Wo{t), 
Zo{t) are present. That is because, at j ^ ji-es, use of the recursion relations H2.13flfl 
and the definition Ij2.f 3d|) of the functions Fj , Gj , Hj , Kj , leads always to expansion 
coefficients Uj, Vj, wj, Zj that are expressed by the zeroth order functions uoit), voit), 
Wo{t), zo{t). 

In the following we shall present the results of the calculation belonging to each 
compatibility condition. For the resonance j — S the calculation can be performed 
easily by hand, but for j — A the computer program Maple had to be invoked in 
order to perform the analytic manipulations. As a result of the Maple program, all 
the coefficients which multiply the higher order powers of (fit proved to be analytically 
zero. The expression associated with the zeroth order power of (j)t has been evaluated 
further by hand to get the final results which will be presented and discussed below. 



3.1. Gonditions for the potentials arising from j = 3 

In section 12 it has been shown that at j = 3 the compatibility condition decomposes 
into two distinct parts which are however not independent of each other [see equations 
H2.20qI -6'1 and the remark thereafter]. 

The elaboration of the compatibility conditions H2.20al -&') related with j = 3 yields 
the relations 

Fawo - G3M0 = — > - V^i*i)uoi'o + V'lo.oi'o - V\o,oWo = , (3. fa) 

i/szo - A'awo = — > {V2.1 - V2*i)tuo2o + ^"20,0^0 - = . (3.I&) 

It is clear from equation l|2.f 2|1 that only the products Uq^^o ^^^d wqZq are determined 
uniquely so that one element of each pair can be chosen arbitrarily. With the choices 
ito = 1, Wo = I the above relations can only be satisfied if 

^i,i-Ki=0 and 1/10,0 = 1^10,0 = 0, (3.2a) 

^2,1 - = and 1/20,0 = ■^^20,0 = 0- (3.26) 

These conditions show that the expansion coefficients Vi,i and V2,i are real. Moreover, 
using the definition l|2.7() for the expansion coefficients V^.j and the results H3.2afl . we 
get 



f 90^10 (x,t) 

U — l'10,0 — 777 



0! dx^ 



-yio(e(i),t). (3.3) 

=m 



Since this equality holds for any arbitrary function ^(t), it follows that l/io(a;,t) must 
vanish. Similar argumentation leads to disappearance of V2o(a;, t). 

In summary we conclude that in order that the equations (|2.3al -M pass the P-test 
the inhomogeneity terms must vanish and the first order expansion coefhcient of the 
external potentials should be real: 

Fio = 1/20 = 0, (3.4a) 
Vii = Vl^ and 1/2,1 = 1/2*1 • (3-46) 

3. 2. Gonditions for the potentials arising from j = 4 

Without presenting the details of the algebraic manipulations, we state that the 
compatibility condition (|2.22|) (partly with the aid of the formula manipulation 
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^UoWo(Vl,o - Ko)^ + '^wozoiV2M ~ V^o)^ 

-Ti2WoZo{V2,2 + ^2*2) - 72iUoWo(^4.2 + V\*2) 

T21 d T12 d 

-i-^"oWo^(Vi,o - Ko) - '-^i^wozo — {V2fi - V2fl) = (3.5) 

in which, as expected, only the zeroth order coefRcient functions UQ{t)^ voit): Wo(t), 
Z(j{t) appear together with the parameters A, •& and T^j (i,j = 1,2) in a special 
combination. This condition looks much more complicated than that obtained above 
[cf. with equations H3.1al -6)]. Moreover both potentials Vi and V2 are occurring within 
a single relation. 

Let us now write the external potentials in the form 

Vi{x,t) — a{x,t) + ip{x,t), (3.6a) 

V2ix,t) ^-f{x,t)+i5{x,t) {3.6b) 

where a, 7 and /?, 6 are real functions. Exploiting the reality of and V2,i expressed 
by relation l)3.4&|l and using the definition (|2.7|l . we get the results 

(3{x,t)= I3{t) and 5{x,t)=5{t). (3.7) 

This condition which can be checked easily by direct substitution tells that the 
imaginary part of the potential may depend only on the time t but not on the space x 
variables. Using this last result and inserting the definitions H3.6al -&) into the relation 
H3.5(l we get the following expression 

--r^2lU0W0P - ^^12^02^0^' + T^2lW0«0-7r + ^Ti2WqZo — 

A v A at v at 

-Ti2Wozo-g^ -T2iuovo-^ = . (3.8) 

Because the quantity in the square bracket depends only on time t, integration by x 
twice yields the following results 

Ti2Wo-zo7 + T2iUoVoa = Ci{t)x^ + C2{t)x + C^{t) (3.9) 

where the coefficients Ci(i), C2(t), C-i{t) depend only on time and C2, C3 are 
arbitrary real functions. By re-substituting this latter equation into expression (|3.8|) . 
we get the constraint for the function Ci {t) as follows 

T21 /ld/3 .2V ^12 /ld<5 2\ .o.n^ 
C,{t) = —novo [^^-f3 )+ —^ozo (2 d^ " ^ j • (3-10) 

We emphasize that the above result does not mean a restriction for the individual 
form of the real part of the external potentials Vi and V2. As equation (|3.9() shows 
only a weighted sum of the real parts a and 7 is constrained by the compatibility 
conditions (|3.5|l belonging to the resonance j — 4. 

Let us now exhibit a possible consequence of the general constraints (|3.9|l and 
H3.10|l for the potentials by starting from an obvious splitting of the coefficient Ci{t) 
as follows 

Ci{t) = ci'\t)+C['\t) (3.11) 

with 

^ _ ,A . ^ _ ,A ,3,,, 
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Then the compatibility condition expressed by equation H3.9|l can be satisfied by the 
following choices: 

T12W0207 = C[^^x^ + Ci^^x + ^ + f{x,t) (3.13a) 
T2iUQVoa = Ci^'^x^ + C^^'^x + cP ~ f{x,t) (3.13&) 

where c'^^\ C^'^'' {i — 2, 3) are arbitrary real functions of time t and / is an arbitrary 
real function of x and t. Using the expressions 13.6al -6) wc get the following possible 
form for the external potentials: 

where V^''^'', v}^\ V.^"^^ and Vj*'^'' are arbitrary real functions of time t and V{x,t) 
represents an arbitrary real potential function which may be used conveniently in 
design of experiments with BEC. At this point wc have to note that these formulae 
cannot be used for the uncoupled case, since when T12 — T21 = 0, then the 
compatibility condition (|2.22l) changes to (|2.24q1 -6). and in this way V{x,t) does not 
arise. 

In summary we conclude that in order that the coupled GP equations \2.'i(A h) 
pass the P-test, a special combination of the real parts of the potentials Vi and V2 may 
depend only quadratically and/or linearly on the spatial coordinate x. A stringent 
relationship can be established between the coefficient of the quadratic terms and 
the imaginary parts which, in turn, may depend only on time t. An additional real 
potential V of general form may be introduced which explicitly exhibits coupling 
between the external potentials Vi and V2. 



4. Discussion of the results 



In this section we discuss the results from various points of view and make comparison 
with related results obtained by others. 



^.1. Presence of source terms 

In the course of the theoretical study of two-component BEC with attractive in- 
teraction, it has been found |22I that the decay and growth of number of atoms is 
best accounted for by introducing an imaginary contact interaction term in the GP 
equations. We now see that our analysis enables the existence of such source terms, 
by appropriately chosen f3{t) and 6{t) [see equation ()3.14at 6')]. This result holds also 
in the case of one-component BEC as found by Clarkson . 



4-2. Uncoupled case 

Next, we investigate the case T12 — T21 = 0, when the system of the GP equations 
^l.'6(A -h) is decoupled. As an example we derive the resonances. Our general equations 
should reduce twice to earlier results obtained by Clarkson "Si? for the one-component 
GP equation. Starting from the general expression H2.14|l and applying the useful 
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formula (|2.12|) one obtains 

det (Qij)) = X^^fij + 1) fij - mj - 4) - 3j + 4 - 2Zi^^M±^) 

= X'd'U + 1) - - 4) {f - 3j - 4) = X'd'^ij + iffU - if {3 - Af . (4.1) 

The resonances (—1, 0, 3, 4) are those found by Clarkson ^30" and all have a multiplicity 
of two as a result of the double number of the (uncoupled) GP equations. 



4-3. Sign of the potential 

One of the experimental situations where the coupled GP equations (|2.3al -&') serve as 
a theoretical basis is the creation of two component BEG |16| . In such experiments 
alkali atoms are confined by symmetrically arranged harmonic trap potentials. One 
of the possibility of our results is to choose in equations H3.14qI -6') all functions Vi^\ 
V^'^\ ^2^"'^', V^j'"', (3, S equal to zero and let the potential V{x,t) operate as a field 
confining the alkali gas particles. It is then required that in equations H3.14fll -&) the 
terms in which our confining potential V occurs do have the same sign. The condition 
that those two terms with V have the same sign is in general 

T12T21 (AT22 - m2) {mi - AT21) < (4.2a) 

which can be expressed also in terms of the scattering lengths as 

012021 (Aa22 - i?ai2) {dan - Xa2i) < . (4-2&) 

Because, physically 012 — 021, the above condition for the equality of signs of the V 
terms in equations H3.14qI -&') is fulfilled for the usual experimental case with X — d ii 

On < 012 < 022 or a22 < 012 < on. (4.3) 

If the scattering lengths ai2 = 021 are greater or lesser than both an and 022 then the 
sign of the terms containing the arbitrary potential V is different, which corresponds 
to untrapping one of the BEG components. 



4- 4- Fulfillment of equation 1^2. 18\) 

The best studied example of the two-component BEG involves Rb atoms in two 
different hyperfine states. It has been found experimentally [71IH1EI) and numerically 
[T^ that a stable configuration of soliton-like vortex in the two-component condensate 
is achieved in the case where the scattering lengths are in the proportion: 

an : 012 : 022 — 1-03 : 1 : 0.97, with ai2 = 021 . (4-4) 

Let us now check whether these ratios obey our general condition (|2.18|l with integer 
m. Since X/d = 1 expression (|2.18() can be written as follows 

2(oii/a2i)(o22/ai2) - (011/021) - (022/012) _ J_ 

(0ll/a2l)(022 /0i2) - 1 16 

The insertion of the above ratios gives 

2 -1.03 -0.97 - 1.03 - 0.97 -0.0018 1 ., n , ^ 

= = 2 EE— 2m +1)2 + 7 (4.6) 

1.03-0.97- 1 -0.0009 16 ^^ ' ^ ^ ' 

which yields 

m = 2. (4.7) 



[(2m +1)2 + 7]. (4.5) 
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This result means that the experimental ratios (|4.4|) correspond to just one of the 
possible solutions of the GP equations characterized by a m = 2 potential family. 

Proceeding further, one can determine the resonances belonging to the experi- 
mentally found ratios H4.4|l to be [cf. with equations (|2.1t)|) . H2.17|l ] 

ji=2 + m = 4, j2 = l-m = -l. (4.8) 

This result means that no further work is needed, the underlying potential 
falls into the category defined by the compatibility condition for j — 4; a, possible 
representation of such potentials is given by t\'6.l4:d -b). Indeed, the quadratic trap 
potential used in the experiments suits well to the general form of potentials obtained 
from the analysis of the resonance at j = 4. 

5. Summary 

In this paper the first step towards verification of integrability of the coupled GP 
equations by means of the P-analysis has been carried out. It has been shown 
that the GP equations pass the P-test provided a special relation among the system 
parameters (masses, interaction strengths) is satisfied [cf. with (|2.18ll ]. One of the 
recent experiments has been taken as an example. In this experiment [71 |H1 El and 
a subsequent numerical study ^| , the vortex stability of a two component BEG has 
been investigated. It is found that the system parameters at which stability occurs are 
just in the proportion which fits our relation H2.18f) with m = 2, a condition necessary 
for the GP equation to pass the P-test. 

As the GP-equations play a great role in describing BEG, a particular attention 
has been paid to establish the admissible forms of the confining trap potentials of 
experimental interest. It has been found that, in addition to the prescribed form 
resulted by the P-analysis of a single GP equation, there is a possibility of introducing 
an extra potential term of arbitrary shape into the external potentials [cf. with 
equations (|3.14qI -6)]. Also, some discussion of the results has been added which 
includes the comparison of the earlier results obtained for the one component GP 
equation, the role of the source (imaginary) terms /3(t) and 6{t) in the potentials, and 
the sign of the additional potential terms. 

Finally, we add a remark to the fulfillment of equation H2.18f) with integer m's. 
In the light of experimental errors the above agreement m — 2 may seem to be 
accidental. We note however that soliton-type structures (e.g. vortices in 3D) possess 
an outstanding stability sometimes called 'robustness' which enables these particle- 
like formations to survive for a long time or even to arise in circumstances that do 
not fit the exact constraint of mathematics. Therefore equation H2.18|l may prove also 
useful in exploring other regions of parameters where such stable structures are to be 
observed in binary condensates. 
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